1/1 


DVNAAICAL  ANALYSIS  OF  NOLECULAR  DECAY  AT  SPHERICAL 
SURFACES(U)  STATE  UNIV  OF  NEH  YORK  AT  BUFFALO  AHHERST 


UNCLASSIFIED 


P  T  LEUNG  ET  AL. 


AUG  87  TR-38  N88814-8«-K-8i4J 

F/G  7/4 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARD*  1963-A 


AD-A184  908 


OTIC  FILE  COEJ 


OFFICE  OF  NAVAL  RESEARCH 
Contract  N0001 4-86-0043 
TECHNICAL  REPORT  No.  50 


Oynamlcal  Analysis  of  Molecular  Decay  at  Spherical  Surfaces 
by 


P.  T.  Leung  and  Thomas  F.  George 

Prepared  for  Publication 
In 

Journal  of  Chemical  Physics 

Departments  of  Chemistry  and  Physics 
State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


DTIC 

electee 

SEP  2  3 19871 


August  1987 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 
United  States  Government. 

This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


y. 


I 

V- 

iV 

# 

r 
ii 
$ 
; i 

v 


$ 

& 

k 

b 


7  09  15  ...  1®  5 


i 

I 

I 

1 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE 


1*.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2A  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


2b.  OECLASSIFICAT I ON/OOWNG  RACING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERIS) 


UBUFFAL0/DC/87/TR-50 


6a  NAME  OF  PERFORMING  ORGANIZATION 

Depts.  Chemistry  4  Physics 
State  University  of  New  York 


6c.  AOORESS  (City.  Stale  and  /.If  Coda) 

Fronczak  Hall,  Amherst  Campus 
Buffalo,  New  York  14260 


3.  DISTRISUTION/AVAILASILITY  of  REPORT 

Approved  for  public  release;  distribution 

unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBERIS) 


b.  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
lit  applicable) 


7b.  AOORESS  icily.  State  and  ZIP  Coda ) 

Chemistry  Program 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


Ba  NAME  OF  FUNOINGISPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
II f  applicable) 

Contract  N00014-86-K-0043 


8c.  AOORESS  (City,  Stole  and  '/.IP  Code ) 

|  tO.  SOURCE  OF  FUNDING  NOS. 

Chemistry  Program 

PROGRAM 

PROJECT  i 

task 

WORK  UNIT 

800  N.  Quincy  Street 

Arlington,  Virginia  22217 

ELEMENT  NO. 

NO. 

NO. 

NO. 

11.  title 

Dynamical  Analysis  of  Molecular  Decay  at  Spherical  Surfaces 


12.  PERSONAL  AUTHORIS) 

P.  T.  Leung  and  Thomas  F.  George 


14.  DATE  OF  REPORT  lYr..  Mo..  Day ) 

August  1987 


16.  supplementary  notation 


13a  type  of  report 


13b.  TIME  COVEREO 
from _  TO. 


15.  PAGE  COUNT 

13 


Prepared  for  publication  in  Journal  of  Chemical  Physics 


COSATI  COOES 


GROUP  SUB.  GR 


IS.  SUBJECT  TERMS  (Continue  on  rrpfrif  i/ n«ct«Mry  and  identify  by  block  number ) 

MOLECULAR  DECAY  DIELECTRIC  SPHERE 

DYNAMICAL  ANALYSIS  DIFFRACTION  THEORY 

SPHERICAL  SURFACES  DIPOLE  ANTENNA 


ttnue  on  reverie  tf  neceeemry  end  identify  by  block  number t 

A  dynamical  approach  to  the  classical  decay  rates  for  molecules  near  a 
dielectric  sphere  is  presented  through  the  application  of  the  diffraction  theory 
for  a  dipole  antenna  established  by  Van  del  Pol  and  Brenner.  This  theory  is 
somewhat  simpler  than  but  formally  equivalent  to  that  established  by  Ruppin  and 
preserves  a  feature  which  is  closer  to  the  method  of  the  theory  established  by 
Chance,  Prock  and  Si  1  bey  for  a  flat  surface.  The  results,  when  compared  to  those 
obtained  from  the  static  image  theory,  show  that  this  latter  theory  can  be  very 
inaccurate  for  large  molecule-sphere  distances  or  highly-conducting  spheres, 
consistent  with  previous  findings  for  surfaces  with  perfect  flatness  or  small 
roughness. 


2a  OlSTRIBUTIONIAVAl LABILITY  OF  ABSTRACT 
UNCLASSIFIEO/UNLIMITEO  S3  SAME  AS  RPT.  8!  OTIC  USERS  □ 


22a  NAME  OF  RESPONSIBLE  INOIVIOUAL 

___  Dr.  David  L.  Nelson 


DO  FORM  1473, 83  APR  eoition  of  i  jaf 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22b.  TELEPHONE  NUMBER  22c  OFFICE  SYMBOL 

< Include  Area  Code) 


(202)  696-4410 


EOITION  OF  1  JAN  73  IS  OBSOLETE. 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


DYNAMICAL  ANALYSIS  OF  MOLECULAR  DECAY  AT  SPHERICAL  SURFACES 


P.  T.  Leung  and  Thomas  F.  George 
Departments  of  Chemistry  and  Physics  &  Astronomy 
239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
Buffalo,  New  York  14260 


ABSTRACT 


A  dynamical  approach  to  the  classical  decay  rates  for  molecules  near  a 
dielectric  sphere  is  presented  through  the  application  of  the  diffraction  theory 
for  a  dipole  antenna  established  by  Van  del  Pol  and  Brenner.  This  theory  is 
somewhat  simpler  than  but  formally  equivalent  to  that  established  by  Ruppin  and 
preserves  a  feature  which  is  closer  to  the  method  of  the  theory  established  by 
Chance,  Prock  and  Silbey  for  a  flat  surface.  The  results,  when  compared  to  those 
obtained  from  the  static  image  theory,  show  that  this  latter  theory  can  be  very 
inaccurate  for  large  molecule -sphere  distances  or  highly-conducting  spheres, 
consistent  with  previous  findings  for  surfaces  with  perfect  flatness  or  small 
roughness . 


I.  INTRODUCTION 


The  theoretical  aspect  of  the  problem  of  molecular  lifetimes  in  the  vicinity 
of  a  flat  or  rough  dielectric  (often  metallic)  surface  has  recently  received 
considerable  attention  (see,  for  example.  Refs.  1*3).  It  is  known  that  for 
molecules  a  few  Angstroms  away  from  the  surface  so  that  the  "quantum  spreads"  of 
both  the  surface  electrons  and  molecule  can  be  neglected,  the  classical 
phenomenological  (CP)  approach  works  adequately.  Within  the  CP  approach,  there 
are  the  exact  dynamical  (energy  transfer)  theory  (ET)  and  the  approximate  but 

simpler  static  image  theory  (IT).  In  spite  of  the  fact  that  the  exact  theories 

(ET)  have  been  available  for  some  time  now  for  both  the  cases  of  flat1  and 

spherical^  surfaces,  the  simpler  theories  (IT)  have  still  been  applied  many  times 

6  7  2  3 

to  cases  of  fla*,  spherical  and  other  kinds  of  rough  *  surfaces. 

Justification  has  then  been  given  based  on  the  argument  that  provided  the 

molecule-surface  distance  (d)  is  much  shorter  than  the  emission  wavelength  (A), 

17  8  9 

IT  should  be  as  accurate  as  ET.  *  *  In  a  recent  paper,  however,  we  clarified 
(with  reference  to  a  flat  surface)  that  the  condition  d  «  A  is  not  sufficient 
and  that  IT  can  be  very  inaccurate  for  highly-conducting  substrates,  even  though 

such  a  condition  is  realized  as  in  most  experimental  situations.  In  a  subsequent 

paper, 1®  we  also  established  a  dynamical  theory  (ET)  for  rough  surfaces. 
Moreover,  the  theory  established  in  this  latter  paper  is  a  perturbative  theory 
which  is  restricted  to  be  practical  only  for  very  shallow  roughness. 111 

There  remains,  therefore,  the  problem  of  a  dynamical  theory  for  molecules 
decay  on  a  surface  of  large  roughness.  For  this  case,  however,  it  has  been  found 
that  a  tractable  model  is  obtained  by  replacing  the  rough  surface  by  a  collection 
of  spheres  (or  spheroids)  and  allowing  the  radii  of  the  spheres  to  be  arbitrarily 
large.  Usually,  in  this  approach,  the  exact  solution  for  an  isolated  sphere  is 
worked  out  and  then  a  cluster  of  neighboring  spheres  is  considered  to  model  the 


I 
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actual  surface.^’11  Nevertheless,  in  aost  of  this  previous  work^*11,1^  except 
that  by  Ruppin,^  IT  has  been  applied  with  again  the  erroneous  justification  for 
conditions  satisfying  d  «  A. 

In  this  paper,  we  shall  reformulate  the  problem  for  the  dynamical  decay 
rates  of  an  admolecule  near  a  dielectric  sphere  in  a  different  but  simpler 
approach.  Our  main  goal  here  is  to  compare  the  dynamical  treatment  with  the 
static  one  (IT)  and  to  point  out  that  this  latter  thoery  can  be  very  inaccurate 
for  highly-conducting  spheres  in  spite  of  d  «  A.  Since  the  formal  theory  (in  a 
different  approach)  has  already  been  available,^  here  we  shall  present  the  detail 
of  only  the  simplest  case,  namely,  a  perpendicular  dipole,  for  the  sake  of 
illustrating  the  points  we  have  addressed  above. 

II.  VAN  DEL  POL-BREHMER  THEORY 

In  spite  of  the  theory  worked  out  by  Ruppin5  who  applied  the  dipole 

scattering  theory  of  Kerker  et  al1^  which  is  in  turn  based  on  the  Lorenz-Mle 
14 

theory,  one  always  wonders  whether  such  a  problem  can  be  formulated  in  a 
fashion  which  is  closer  to  the  original  dynamical  theory  established  by  Chance, 
Frock  and  Silbey  (CPS)  for  flat  surfaces1  through  the  application  of  the 
Soesnerfeld  theory  for  radiating  dipoles  above  a  "flat  earth". ^  Indeed,  some 
time  after  Sommerfeld  published  his  work.  Van  del  Pol  and  Brenner  had  generalized 
Sooner f eld' s  problem  to  the  case  of  a  spherical  earth  and  had  shorn  that  the 
Sommerfeld  theory  is  recovered  in  the  limit  where  the  radius  of  the  earth  becomes 
infinitely  large. To  apply  their  theory  to  our  problem,  we  recall  that  the 
Berts  vector  for  the  region  outside  the  sphere  is  given  (in  spherical 
coordinates)  by1 ^ 
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where  we  have  assumed  a  molecular  dipola  p  located  at  (d,0,0)  and  oriented 

radially  above  a  aphere  of  radiua  a.  k  is  the  emission  wave  number,  R  is  the 

distance  measured  from  p,  and  J  .  h^  and  f  are  the  usual  spherical  Bessel 

n  n  n 

functions  and  Legendre  polynomials,  respectively.  The  spherical  reflectance  Rr 


in  Eq.  (1)  is  given  by 


16 


(l+n)(l-e)  +  eka[Jn+1(ka)/Jn(ka)  -  Jn+1(/eka)//tJn(/cka)] 
n  "  (l+n)(c  1)  ♦  eka[jn4l(/Eka)//cjn(/eka)  -  hj^kaj/h^tta)] 


(2) 


where  c(u)  is  the  complex  dielectric  constant  of  the  substrate  sphere.  By 
considering  only  the  reflected  field  and  using  the  expression  for  the  radial 
electric  field*® 


B  (r.8)  -  (k2  ♦  ~)(f  D)  , 
r  ir2  d 


(3) 


we  obtain  finally  the  reflected  field  at  the  dipole  site  in  the  form 


I  <d,0)  -  ipk  V  R  -M—  lh^1)(kd)]2  n(n+lK2n+l2  # 

nto  hiU<k*>  d 


(4) 


and  tha  G- function  in  this  dynamical  theory  as 


E  (d,0) 

C-_(«)  -  — -  .  (5) 

ET  \t 

Tha  total  (i.a.,  both  nonradiativa  and  radlativa  transfer)  molecular  decay  rata 

1-3 

in  tha  presenca  of  tha  sphara  is  then  given  by 

'kt  *  V1  +  I  jji  Ih®et)  •  (6) 

with  yq  being  the  rata  for  a  free  molecule  and  q  the  quantum  yield  of  the 
emitting  state,  respectively.  Tha  results  obtained  in  Eqs.  (4)-(6)  are 
mathematically  equivalent  to  those  obtained  by  summing  Eqs.  (27)  and  (33)  in 
Ruppin's  paper,  except  that  Eq.  (4)  here  is  derived  and  expressed  in  a  somewhat 
simpler  manner,  since  it  does  not  contain  any  integrals  involving  the  Bessel 
functions  and  there  is  only  one  complex  reflectance  coefficient  appearing  in  the 
final  expression.  We  have  checked  numerically  that  our  results  have  reproduced 
identically  the  results  in  Fig.  1  of  Ruppin's  paper.  For  the  case  of  parallel 
dipoles,  the  results  can  be  obtained  similarly  by  introducing  two  Hertz  vectors^ 
and  again,  one  expects  somewhat  simpler  results  obtained  as  compared  to  Ruppin's 
Eqs.  (28)  and  (34).  In  the  following,  we  shall  compare  Eqs.  (4)-(6)  with  the 
static  image  theory  and  assess  better  the  limiting  case  provided  by  this  latter 
theory. 

III.  LIMIT  OF  THE  IMAGE  THEORY 

It  has  been  widely  argued  that  for  d  «  X,  E?  in  Eq.  (4)  can  be  replaced  by 
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the  static  image  field  given  by 


6 


E;(d,0) 
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where  tha  n-pola  polarizability  an(u)  is  given  by 
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The  G-function  may  then  be  defined  accordingly  (G^  *  EVp),  and  YjT  will  be  just 

as  in  Eq.  (6)  with  G^  replaced  by  G^.  We  have  carried  out  numerical 

calculations  for  both  the  distance  dependence  and  the  frequency  spectrum  of  both 

IfflGgj,  and  ImG^  for  a  sphere  of  radius  a  -  100  A.  Figure  1  shows  the  distance 

variation  of  ImG  according  to  both  ET  and  IT  at  u  *  2.5eV  (X  -  5000  A)  for  both  a 

19 

silver  and  a  nickel  sphere.  It  is  not  difficult  to  see  that  under  these 

conditions  where  d  and  a  are  much  smaller  than  X,  IT  can  be  very  inaccurate  for  a 

highly-conducting  sphere  such  as  Ag,  though  for  the  case  of  a  Ni  sphere,  IT  and 

ET  are  fairly  close  to  each  other  for  this  range  of  distances.  This  is 

9  10 

consistent  with  the  previous  observations  for  a  flat  and  a  shallow  grating 

surface,  and  the  physical  origin  for  such  a  phenomenon  has  been  well  explained  in 

the  previous  papers.  *  Figure  2  shows  the  frequency  spectrum  of  ImG  for  d  ■ 

500  A  for  a  Ag  sphere,  for  which  IT  is  expected  to  break  down  appreciably. 

Nevertheless,  for  the  small  region  close  to  the  surface  plasmon  resonance  (-3.5 

eV),  the  relative  agreement  between  ET  and  IT  is  the  best.  This  is  in  contrast 

to  the  previous  comparison  for  the  shallow  grating  case*®  and  may  be  due  to  the 

20 

fact  that  for  ka  «  1,  both  ET  and  IT  have  a  very  similar  resonance  structure. 


IV.  CONCLUSION 


In  this  paper,  we  have  presented  an  alternative  formulation  of  the  dynamical 

decay  rates  of  molecules  near  a  spherical  surface  following  an  approach  which 

differs  from  that  of  Ruppin^  and  stays  closer  to  the  framework  of  the  original 

CPS  theory  for  a  flat  surface.^  In  fact,  it  is  straightforward  to  show  that  Eqs. 

(4) -(6)  lead  back  to  the  results  for  a  perpendicular  dipole  in  the  CPS  theory  by 

taking  the  limit  a  -*■  ■>,  d  •*  •  ,  but  with  the  difference  d-a  kept  as  a  finite 

constant.*®  Furthermore,  in  spite  of  the  availability  of  the  Van  del  Pol-Bremmer 

theory  which  takes  the  advantage  of  the  concept  of  the  Hertz  vector,*®  it  is 

5  13  20  21 

interesting  to  note  that  most  of  the  previous  dynamical  theories  *  *  ’  for 

14 

molecule -sphere  interactions  are  lased  on  the  Lorenz-Mie  theory  whose 
mathematical  structure  is  in  general  more  complciated.  In  the  light  of  the 
present  investigation,  it  seems  that  an  alternative  approach  to  all  these 
previous  problems  based  on  the  Van  del  Pol-Bremmer  theory  is  worthwhile  because 
of  its  comparatively  simpler  structure,  as  illustrated  by  the  sample  calculation 
in  this  paper.  In  addition,  due  to  the  fact  that  the  present  approach  is  more  of 
a  scalar-type  expansion  (in  contrast  to  the  use  of  the  vector  harmonics  in  the 
other  theories),  one  may  find  it  easier  to  generalize  the  theory  to  the  case  of  a 
cluster  of  spheres. 

We  have  further  compared  this  theory  with  the  static  theory  and  have  shown 
once  again^’*®  that,  in  contrast  to  many  previous  expectations, ^'®* ** *  ^  the 
static  theory  can  be  very  inaccurate  for  highly-conducting  (e.g.,  Ag)  spheres,  in 
spite  of  the  fact  that  d  «  A.  Hence,  all  the  previous  work  on  SERS, 
fluorescence  and  other  resonant  absorption  processes^  which  has  utilized  image 
fields  in  their  formalisms  becomes  inaccurate  subject  to  the  present 
observations,  and  therefore  must  be  reformulated  by  introducing  a  dynamical 
description  for  the  decay  rates  of  the  admolecules. 
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FI CURE  CAPTIONS 


1.  Comparison  between  the  energy  transfer  theory  (FT,  solid  curves)  and  the 
image  theory  (IT,  dotted  curves)  for  a  (a)  Ag  and  (b)  Nl  sphere  at  <*  •  2 . 5  eV 
for  a  range  of  molecule-sphere  distances.  The  unit  of  C  Is  A  . 

2.  Comparison  between  ET  and  IT  for  the  frequency  spectrum  ImG(u)  at  d  *  $00  A 
for  a  Ag  sphere. 
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